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ABSTRACT: Interrelation of the Coleman's representabilty theory for 1- 
density operators and abstract algebraic form of the Hohenberg-Kohn theorem is 
studied in detail. Convenient realization of the Hohenberg-Kohn set of classes of 
1-electron operators and the Coleman's set of ensemble representable 1-density 
operators is presented. Dependence of the Hohenberg-Kohn class structure on the 
boundary properties of the ground state 1-density operator is established and is 
illustrated on concrete simple examples. Algorithm of restoration of many electron 
determinant ensembles from a given 1-density diagonal is described. Complete 
description of the combinatorial structure of Coleman's polyhedrons is obtained. 
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I. Introduction 

In the middle of the twentieth century there appeared two papers that 
became contemporary classics and called into being a new branch of quan- 
tum chemistry. The first paper, dedicated, at first glance, to solution of a 
very abstract problem of representability of 1-density operators by ensem- 
bles of many electron states, was written by Coleman and published in 1963 
PP. One year later, in 1964, Hohenberg and Kohn published their theo- 
rem j2] about correspondence between representable densities and external 
potentials. From mathematical point of view both results belong to the func- 
tional analysis and have close interrelation. Coleman's theorem gives explicit 
analytic description of the set of all ensemble representable 1-density opera- 
tors and Hohenberg-Kohn theorem suggests implicit way of parametrization 
of representable 1-density functions by certain classes of 1-electron operators 
(potentials). Further development of Coleman's approach was mainly concen- 
trated on attempts to get the necessary and sufficient conditions of ensemble 
representability for density operators of higher order, especially for 2-electron 
density operators (see [3 EI and references therein). It turned out to be a 
very complicated mathematical problem and to the best of our knowledge 
the constructive analytic description of the set of all representable 2-density 
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operators is still not found. On the contrary, mathematically almost triv- 
ial Hohenberg-Kohn statement became a banner of the density functional 
community and is considered at present as the foundation of the density 
functional theory (DFT) (see, e.g., [S]-[E])- 

The present paper is an attempt to give reasonably rigorous analysis 
of both the Coleman's theory and abstract form of the Hohenberg-Kohn 
theorem for the case when one-electron sector of the Fock space is of finite 
dimension. 

Section II is dedicated to 1-density operators and their properties. Be- 
sides formulation of different versions of the Coleman's theorem, a useful 
realization of the set of all ensemble representable 1-density operators is 
given. Combinatorial structure of Coleman's polyhedrons is discussed in Ap- 
pendix A. In Appendix B an example of constructive description of the set 
of 1-density operators representable by pure states may be found. 

In Section III general iteration formula for restoration of many electron 
determinant ensembles from some fixed representable 1-electron diagonal is 
discussed. This formula first appeared in our electronic publication [§]. Con- 
crete numerical examples of using this formula are given. 

In Section IV a very general algebraic formulation of the Hohenberg- 
Kohn statement is analyzed. Abstract considerations are accompanying by 
concrete examples demonstrating dependence of 1-density diagonal matrix 
elements on parametrized classes of 1-electron operators. 

II. General Properties of 1-Density Operators 

In coordinate representation density operator of the first order (1-density 
operator) is introduced as the integral operator 



acting on the one-electron sector of the Fock space. Here x is a vector of 
space-spin variables and p(x,x') is a kernel of the operator (JTJ). This oper- 
ator is said to be representable by a pure p— electron state *&(xi, . . . , x p ) 
if its kernel can be obtained by contraction (integration) of the product 
. . . , x p )^*(x', . . . , x p ) over p — 1 variables: 



p(x, x') — / . . . / ^(x, x 2 , . . . ,x p )^*(x', x 2 , . . . ,x p )dx 2 . . . dx p (2) 




(1) 
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Contraction of an ensemble of arbitrary finite family of p— electron states 




p(x,x') = ^ A; 



, x p )ty*(x', x 2) ■ ■ ■ , x p )dx 2 . . . dx p (3) 



i=l 
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£A j = 1 
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gives a kernel of the so-called ensemble representable density operator. The 
diagonal part p(x) = p(x, x) of this kernel is called (ensemble) representable 
density function and plays fundamental role in the density functional theory. 

Directly from definition three properties of representable 1-density oper- 
ators are easily deduced: 

(1) = p (Hermiteancy); 

(2) p > (positive semi-definiteness) ; 



Density operators satisfying conditions (l)-(3) are not necessarily rep- 
resentable. The necessary and sufficient conditions of the ensemble repre- 
sentability of 1-density operators was found by Coleman PJ 0] and may be 
formulated as follows. 

Theorem 1 (Coleman). 1-density operator p is representable by an ensem- 
ble of p— electron states if and only if 



Here I is the 1-electron identity operator. 

Now let us turn to the algebraic version of the representability problem 
confining ourselves to the finite-dimensional case. We suppose that 1-electron 
sector of the Fock space Fn,i is spanned by n = \N\ orthonormal molecular 
spin-orbitals (MSOs) {ipi}ieN arid that standard creation-annihilation opera- 
tors are associated with this MSO basis set. The electronic space correspond- 
ing to p— electron system is just the pth Grassmann power Tn, p = [\? Fn,i 
of the 1-electron space. Following Dirac, we identify the algebra of linear 



(3) Tr p = 1. 



(1) 
(2) 

(3) 
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operators over p— electron space J~n, p with the tensor product J~n, p <8> J~n p 
of this space and the space of |-linear functional on Tn, p - The contraction 
operators is defined as 

n 

c — y^aj ® a\. (5) 

i=l 

and it is easy to show that 1-density operator representable by a pure 
p— electron state |^) is 

p=-c p - 1 \^)^\ (6) 
pi 

Due to linearity of the contraction operator Eq.fJUJ) is easily generalized to 
treat ensembles of p-electron states. The set of all representable by ensembles 
of p-electron states 1-density operators will be denoted by the symbol Vn, p ,i 
and will be referred to as Coleman's set. It is a compact (and, consequently, 
closed) convex subset of the affine hyperplane Tr p = 1 situated in the 
operator space Fn,i ® J~nv From Coleman's theorem it readily follows that 
the set Vjv, P ,i is invariant with respect to transformations u <S> where u G 
U(JF iV1 ) is 1-electron unitary operator (the so-called unitary invariance of 
the representability problem). In other words, Vn, p ,i is a U(JF /V1 ) -space with 
respect of 1-electron unitary group action (see, e.g. [T31 PH] ). 

The structure of the set JV, P ,i of 1-density operators representable by 
pure p-electron states is more complicated. It is the image of the projective 
space P(jF/vp) over the p-electron sector of the Fock space with respect to 
the contraction ^c p ~ 1 . Since P(jF A r p ) is compact and connected space, its 
image with respect to continuous mapping (contraction) is also compact and 
connected. In addition, projective spaces do not admit global parametriza- 
tions being the simplest examples of the so-called analytic manifolds (see, 
e -g- ; [U1IJ2])- It may therefore be expected that global parametrization does 
not exist also for the set JV, p ,i as a whole. 

In non-relativistic quantum chemistry molecular spin-orbitals are pre- 
sented as (tensor) products of spatial and spin functions and it is normally 
assumed that p-electron states under consideration have a fixed projection 
Ms of the total spin. 1-electron sector of the Fock space can be decomposed 
into a direct sum of its a- and /5-subspaces and 1-density operator becomes 
a direct sum of a- and /3-components. There are two equivalent forms of 
the representability conditions in this case. If the total 1-density operator is 
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written as p = p a © pp then p is representable if and only if 
(!) Pi = Pa, 

(2) 0< P<T <±I ff ; (7) 

(3) Tr p CT = p f. 

If p is written as p = ^-p a © then it is representable if and only if 

(1) Pt = Pa] 

(2) 0<p a <±I a , (8) 

(3) Tr Pa = 1. 

Here I a is the identity operator in cr-subspace of 1-electron space, p a is the 
number of cr-electrons, and a = a, (3. In the second case p a G Vm,pb,i where 
M is the orbital index set. 

Each 1-density operator is Hermitian and can therefore be diagonalized 
and its eigenf unctions (the so-called natural MSOs) constitute a basis of 1- 
electron sector of the Fock space and without loss of generality this basis can 
be considered to be orthonormal. In Dirac's notations spectral resolution of 
1-density operator p is 

p = ^A l |^)(^| (9) 

where Aj are the so-called natural occupancies. For 1-density operators of 
the form of Eq.fJH} Coleman's theorem can be formulated as 

Theorem 2 (Coleman). 1-density operator p presented as its spectral res- 
olution is representable by an ensemble of p- electron states if and only if 

(1) \i E R for all ieN 

(2) < A, < \ for all % e N; (10) 

(3) E A* = 1. 

ieN 

In geometric terms for a fixed n-frame ip = (ipi, . . . , ip n ) (orthonormal 
ordered MSO basis) the set V]v,p,i °f representable 1-density operators of the 
form (JHJ) is the intersection of the standard simplex and a cube (with the edge 
length -) being therefore a convex polyhedron (situated in the hyperplane 
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Aj = 1). Both parametric and analytic descriptions of this polyhedron 

ieN 

are available. 

Parametric description: Polyhedron V% pl is a convex hull of (™) vertices 

v pll (R) = -J2\^i)(H (11) 

where i? is p-element subset of the MSO index set N. 

Analytic description: It is given by the Coleman's conditions (l)-(3) of 
Theorem El In geometric terms the polyhedron x has In hyperfaces with 
normals 

n° = p\mH (12) 

n, 1 = -p\^)m + J2\^^\ ( 13 ) 

It is clear that for any fixed n-frame if) the polyhedron V^ pl is homeo- 
morphic to the typical (standard) polyhedron Vjv, p ,i constituted by vectors 
(Ai, . . . , A n ) G W 1 satisfying the Coleman's conditions (l)-(3) of Theorem El 
Combinatorial structure of this polyhedron is outlined in Appendix A. 

In terms of natural occupancies it is easy to describe the border of the set 
Vn,p,i' 1^ i s constituted by 1-density operators that have at least one natural 
occupancy equal to zero or I /p. For the border of Vn,p,i standard symbol 

o 

9Vn,p,i will be used. Interior part of Vn,p,i will be denoted as Vn,p,i- Note 
that the minimal number of non-vanishing natural occupancies is equal to 
the number of electrons p and, if it is the case, then all these occupancies are 
equal to 1/p (single determinant density operators that can be considered as 
vertices of Vn,p,i)- 

Now we can describe two different realizations of the set of ensemble rep- 
resentable 1-density operators. In the first realization n-frame ip is supposed 
to be fixed and any representable density operator is written in the form 

P= J2 PiMi)W>i\ ( 14 ) 

where matrix (p i: ,) satisfies conditions (l)-(3) of theorem 1. The set Vn,p,i is 
the image of the mapping 

r : (A, u) — > u\v) (15) 
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where A G V]^ pl and u is 1-electron unitary transformation. It is easy to 
see that this mapping is not injective. In this realization Vn,p,i is a union 
of orbits with respect to the unitary group action where for each orbit its 
representative, diagonal in the basis ip, is selected. 

In the second realization the set of ensemble representable 1-density op- 
erators is considered as a disjoint union (sum) of fibres x 

where N is the manifold of all orthonormal n-frames of 1-electron space !Fn,i- 
In more formal terms Vn,p,i is a total space of (trivial) fibre bundle with the 
manifold TsT as its base and the polyhedron Vjv, p ,i as its typical fibre (see, e.g., 
[TTl 112] ). In such a realization the set of ensemble representable 1-density 
operators is homeomorphic to the Cartesian product N x Vzv,p,i- 1-electron 
unitary group acts transitively on the base N as 

f = mfj = (^ 1 ,...,^ n )U (17) 

where U is a unitary matrix connecting two n-frames 

n 

^ = U ^ = J2^U JU (18) 

and the induced mapping u <8> v) maps fibre x onto fibre v 

The set of 1-density operators representable by pure p-electron states can 
also be considered disjoint union of fibres 

aW = U p w ( 19 ) 

where for each ip e INT the fibre P% pl is a compact connected subset of x 

n 

necessarily containing all its vertices and the central point - I^XV^I- An 

n i=i 

example of explicit description of a fibre P% p x is given in Appendix B. It is 
pertinent to note that implicitly fibre bundles appeared in quantum chemistry 
at its infancy in disguise of the so-called multi-configuration self-consistent 
field (MCSCF) theory. MCSCF fibre bundles have the set N of n-frames as 
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their base (SCF part) and unit spheres of some fixed dimension as the typical 
fibres (MC part). 

If n-frame ip = (/0i, • • • , tp n ) is fixed then for any integer q = 1, 2, . . . , n it 
is possible to define the diagonal mapping 

(?) (?) 
: A Rs\R)(S\ - A RR\ R )( R \ ( 2 °) 

R,ScN ReN 

where symbol \R) stands for g-electron determinant built on spin-orb itals 
ipi , . . . , ip n with indices from subset R. If R = {ri < . . . < r q } then \R) may 
be interpreted as either a factorable p-vector ip ri A ... A ip r in Grassmann 
algebra (electronic Fock space) or as a vector . . . aj |0) obtained by suc- 
cessive application of the fermion creation operators associated with a given 
MSO basis set. 

Definition 1. Vector A G M. n is called representable if it can be realized as a 
diagonal of 1-density operator with respect to some fixed basis. 

It is easy to see that the diagonal mapping (f271jl commute with the con- 
traction operator. This means, in particular, that the diagonal of any repre- 
sentable 1-density operator is representable by an ensemble of determinant 
states and d$ (Vn,p,i) = V% p v Of course, d$(p) G V^ pl does not imply 
P G Vn,p,i- However, the following statement holds true. 

Proposition 1. d^(p) G V]^ pl implies the existence of representable by a 
pure state density operator p' such that d^{p) = d^{p') 

Proof. If d^(p) G V% p i then, by definition, there exists an ensemble of 
p-electron determinant states 

(p) 

A=^2a r \R)(R\ (21) 

RCN 

such that d^(p) = ^c p_1 A, and any p-electron wave function of the form 

(p) 

(0) = exp(iB R )Al\R) (22) 

RCN 
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corresponds to the pure state |\P (0)) (9)\ such that its contraction gives 
representable 1-density operator p' with d^{p'{6)) = d^(p) for arbitrary phase 
vector 0D 

Corollary. 

= V Ar, p ,i (23) 

The connection between two aforementioned realizations of the set of the 
ensemble representable 1-density operators is established by the following 
simple assertion. 

Proposition 2. (1) p G V Np i if and only if d^(p) G V^ pl for any n-frame 

(2) p G 9Vjv, p ,i if and only if there exists n-frame ip such that d^(p) G 

uv N,p,l> 

(3) p G V^ p l if and only if d^{p) G Vn, p ,i for any n-frame tp. 

Proof. Statement (1) readily follows from the Coleman's theorem. As for 
statement (2), it is sufficient to prove that d^(p) G d\/% p l implies p G 9Vjv,p,i- 
Indeed, if p G dV Np i then in the basis tp of natural MSOs d^{p) G <9V^ pl . 
Let us suppose therefore that d^(p) G <9V^ pl which means that pa = 
or pa = i in the MSO basis under consideration. If (Uji) is the matrix of 
coefficients of transformation from this basis to the basis of 1-density operator 
p natural MSOs, then we have 

n 

Pit = ^2 \ Uji\ 2 \ 

where Aj are natural occupancies. If pa = then from this equality it follows 
that there exists at least one index j such that A-,- = 0. If pa = ^ then this 
equality may be recast as 



= £ w (I - ^) 



3=1 

and, consequently, there exists index j such that - — Xj — 0. 

o 

If d^p(p) G p G Vjv, p ,i for any n-frame ip then among natural occupancies 

'N,p,l 



of p there are no occupancies equal to or i. Assumption that p G Vf 
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and there exist n-frame ip such that d^(p) G <9V^ pl contradicts to already 
proved assertion (2) of this Proposition □ 

Explicit description of the border of the polyhedron Vat, Pi i is given in 

Appendix A. Note that d^(p) G V%pi f° r some n-frame ip does not imply 

o 

Iteration formula for generation of all ensembles of determinant states 
that are contracted into a given diagonal is obtained in the next section. 



III. Restoration of Many Electron Determinant Ensembles from 
Diagonal of One-Electron Density Matrix 

With arbitrary vector A*- -* G V^ pl it is convenient to associate three 
index sets: 

Ind(X^) = {teN: A 4 (0) > 0} (24) 

Jndi(A (0) ) = {teN: \f ] = i} (25) 

Jnrf a (A (0) ) = Ind(\ {0) )\Indi(\ {0) ) (26) 

p 

Indices belonging to the last set will be called active. 

Let us present vector A^ G Vff pl as the convex combination 

AW = pp Ro v pll (R ) + (1 - pp Ra )X {1) (27) 
where vertex v p ii(Rq) is defined by Eq.flTTJ), 



xn = E ~r- v w^' m+ E iztss^x*!- (28) 

ieRo FP iEN\R P 

and require the residual vector A^ to be representable. This requirement 
imposes the following restrictions on the admissible values of parameter p R °: 



< < |, * g Ro 

< r& < J, * e N\Ro 



(29) 



The solution of this system is the interval [0, b R °~^ C K where 



6^° = nun{miii{Aj 0) }, min {- - \f ] }}. (30) 

ie.Ro i£N\Ro p 
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If n Ro G (0, b R °] then we arrive at a non-trivial representation of A^ -* as a 
convex combination of vertex v p ii(Rq) and a certain representable residual 
vector A^ 1 ). 

Proposition 3. Arbitrary vector A*- -* G V^ pl admits presentation in the 
form of Eq. \27\) if and only if 

Jndi(A (0) ) C R C Ind{X {0) ) (31) 

Proof. It is sufficient to show that the condition (J31|) is equivalent to 
existence of non- vanishing boundary parameter b Ro . But this follows directly 
from Ea. (13711) □ 

Thus, with each A^ G Vff t we can associate the set V(X^) of 
p— element subsets satisfying the condition (|3Tj). It is easy to show that 

W(X^)\=( naiXi °!L) (32) 



v p-ni(A(°)) 

p 

where n (A^) = \Ind a {X^)\ and m(A^) = |/ndi (A<°>)|. 

Proposition 4. If parameter fi Ro in Eq. fl£7| ) is taken equal to its boundary 
value b Ro then R & V(\ (1) ). 

Proof. If /i fio = b R ° = - then Indi(X w ) = and, obviously, R £ 

' P p 

V(X^). Let us suppose therefore that < b Ro < -. From Eq. (j3TJ]) it follows 
that there exists index G Ind(X^) such that either i* G Ro and jj Ro = 
6«o = A W or ijf e n \ Rq and ^Ro = b Ro = I _ From Eq.(|2SJ it is 

easy to see that in the first case A^ = and, consequently, G i?o but 
2* G' Ind(X^). In the second case A^ = - which means that i» G Jnrfi(A (1) ) 
but z* G" i? - In both cases R & P(A (1) ) □ 

Proposition 5. If parameter fi Ro in Eq. \2T\j is taken equal to its boundary 
value b Ro and this boundary value is different from - then 

Jnrfi(A {0) ) C Indi{\ {1) ) (33) 

V V 

Proof. Let us show that A^° = ^ implies A^ = K Indeed, if fi R ° = 
6% = \(°) < I and \(°) = I for some ; e # then \« = Iz^J = 1. if 
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on the other hand, fi Ro = b Ro — - — and A^ = - for some i G Rq then 



1 — p/j R ° = pXf and A 



(0) 



I □ 



P A<»> V 

Iterating of Eq. pTjl leads to the following expression 



fe-i 



i=0 



'i-l 

U.0--W 

.3=0 



p^ R 'v pll (Ri 



where 



fi Rt G (0, 



k-l 
i=0 



Pfj, 



Ri 



\ {k) (34) 



(35) 



Indi(X {k) ) C i?i C Ind(X (k) ) 



(36) 



for i = 0,l,..., A; — 1. Note that the residual vector in Eq . ()34p is necessarily 
representable. 

Theorem 3. For any vector A^ G V$ 1 the residual vector in iteration 
formula jff^| ) vanishes after at most n a (X ( - ^) — 1 steps if at each step the 
boundary value of parameter fi Ri is selected. 

Proof. Direct consequence of Propositions 1-3. 

Corollary 1. The set V$ pl is the convex hull of (™) vertices v P ii(R). 
Corollary 2. 

Ensemble 



A(Rq, Ri, . . . , Rk f ) — 2j 



i=0 



i-l 



i=o 



p^R^Ri 



(37) 



of p- electron determinant states generated recurrently on the base of the iter- 
ation formula P^D with boundary values of parameters n Ri includes pairwise 
distinct p- element subsets and 



ic p - 1 A(i? ,i? 1 ,...,i4 / ) = A(°). 
pi 1 



(38) 



Definition 2. For X^ G V^ p 1 any its expansion corresponding to the bound- 
ary values of parameters \i Ri will be called boundary expansions. 
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Let us consider simple example of boundary expansion. For N = 
{1, 2, 3, 4} and p = 2 let us take representable rational vector 

V20'4'4'20; 

and select R = {2,4}. With such a selection the boundary value 6 24 = 
and Eq.(j2U) gives 

A(°) = -v P n(24) + -A« 
10 pm ' 10 

where 

AW = f- - - 

\2' r 14' 

The next admissible subset i?i may be chosen, say, as i?i = {1,2}. In this 
case b 12 = | and 

A« = ^(12) + |a« 

where 

A (2) = Q,0,~ Oj =«p;i(13). 
The final expansion (which is by no means unique) is 

A (0) = ^«pu(24) + ^(13) + ^ Pi i(12). 

Note that if G x is rational then all intermediate vectors X^ in 
Eq. (jMj) are also rational if at each iteration the boundary value of parameter 
fi Ri is selected. 

Iteration formula ()34)1 is valid for arbitrary choice of n Ri E (0, but 
Proposition 2 is not. In other words, the conditions ((HSj) do not necessarily 
filter out p— element subsets selected on the previous steps. This may be used 
to construct expansions different from the boundary ones. And even more, 
it is easy to show that on the base of iteration formula (J3"4^ any pre-defined 
p-electron ensemble of determinant states may be obtained. 

Proposition 6. Let 

(p) 

a=J2\Cr\ 2 \R)(R\- (39) 

RcN 
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be an ensemble of p- electron determinant states and A = ^c p X A is the cor- 
responding 1-density diagonal with components 



1 W 

A, = -V|^| 2 . (40) 

Then, using iteration formula \3$ , it is possible to restore the initial ensem- 
ble A. 

Proof. Let us suppose that p-element subsets in Eq. fJHU]) corresponding 
to non-zero coefficients \C R \ 2 are ordered in some fixed (say, lexical) order. 
From Eq. (j4*U|) it readily follows that 

l -\c„\ 2 < x. 



V 



for any i G R. On the other hand, the equality Ai + - \C R \ 2 = - (that is 
a direct consequence of the normalization condition for A) implies that 

-\c R \ 2 <--x i 

p p 

for any i ^ R. Thus, we can guarantee that fi R ° = ^|(7r | 2 G (0,b Ro ] and 
expansion Eq.(|27[) takes the form 

A = A<°> = IC^W^Ro) + (1 - |C* | 2 )A« 

where A*- 1 - 1 corresponds to the ensemble 

(p) 

A (1) = 1 _, r , 2 E \ C n\ 2 \R)(R\ 



Ro\ 



RCN 



Expanding A^ in the form of Eq.(j2ZJ) with fi Rl = ^jq^ py , we come to the 
equality 

A (0) = [C^W^Ro) + \C Rl \ 2 v pll (Ri) + (1 - |C Ro | 2 - \C Rl \ 2 )\( 2 \ 
etc. On the pre-final step k we have the ensemble 

A^ 1 - ) = |C flfc _ i | 2 |-Rfc-i)(-RA:-i| + \C Rk \ 2 \Rk){Rk\ 

14 



where 



\C 



\n 12 



Rk 



k-2 



\C 



\n 12 

I Wife I 



fc-2 



I 2 i-EI^I 2 

i=0 i=0 

Non-zero components of the corresponding vector A^ -1 ' are 



A 



(fc-i) 




2 + \C 



if i e R k -i\Rk 

if Z G R k \Rk-l 

if z e -Rfc H Rk-i 



It can be easily shown that /Z**- 1 = fr^- 1 = ^|C^ fe _J 2 and, consequently, 
= u pi i(i? fe ) □ Now we may suggest the following scheme for generation 

Table 1: Example of application of iteration formula (}3*4"j) for the case N = 
{l,2,3,4},p = 2, and A^ = (£, |, f, 



Interval /i J 






{1,2} 


(o,|] 


1 

5 


(I 1 A I) 

I45 12' 12' 4/ 


1 


{1,3} 


(o,|] 


1 

6 


/l 1 3 3\ 
V8' 8' 8' 8/ 


2 


{1,4} 


(0,11 


1 

16 


/ 1 13 5 \ 
V 14 ' 7' 7' 14/ 


3 


{2,3} 


(o,f] 


1 
7 


VlO' ' 5' 2/ 


4 


{2,4} 


{0} 







5 


{3,4} 


(0,§1 


2 

5 


(1,0,0, |) 



of non-boundary expansions. First it is necessary to arrange p— element 
subsets from Ind(M°>) in some fixed order R , R 1 , . . . , R kf and step by step 
split vertices first v p ±x(Rq) from \(°>, then t> p | 1 (i? 1 ) from A^, etc, selecting 
at each step a non-boundary value of parameter n Ri . This procedure is 
continued till the moment when on some step k the number of remaining 
subsets Rk+i, • • • , Rk f becomes equal to the number of active indices in the 
residual vector A^ minus 1. Then to this residual vector the algorithm 
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of generation of the boundary expansion is applied. Simple example of an 
expansion, obtained in such a way, is given in Table Q for the case 

(20' I' 4' ^) alld P = 2 - 

IV. Algebraic Version of The Hohenberg-Kohn Theorem 

The set of Hermit ian linear operators over p-electron sector J-^ of the 
Fock space will be denoted as "K p . 
The mapping 

v 

X : v -> y*J ® . . . ® J (8w (8) J ... ® J (41) 

defines a basis independent linear representation of the set of 1-electron op- 
erators "Hi by p— electron operators from "K p . Operators X(v) are called 
'1-electron operators acting on p-electron space'. More simple but basis de- 
pendent definition of X(v) may be given in terms of the standard creation- 
annihilation operators associated with some orthonormal MSO basis set. 
Let H e "Hp be some fixed operator over p— electron space and let 

P**( v ) = \^i\v) >< *i\v)\ (42) 
i=i 

be the projector on the eigenspace W^(v) of operator H+X(v) corresponding 
to its lowest eigenvalue E^{v). Let us introduce the following equivalence 
relation on the set 

R H (v,v')^P l H (v)=P 1 H (v') (43) 

Quotient of "K\ modulo the equivalence relation contains classes of equiv- 
alent 1-electron operators and the notation Qh = "Ki/Rh for the set of such 
classes will be used. Class represented by some 1-electron operator v is de- 
noted by the standard symbol [v]. The set of scalar operators is a subset of 
the zero class [0] not necessarily coinciding with it. 

If otherwise is not stated, the quotient topology is supposed to be fixed 
on the set Qh- It is to be emphasize that the definition of the set Qh in the 
case of infinite dimension should be modified to exclude 1-electron operators 
v such that the discrete part of spectrum of the operators H + X(v) is empty. 
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Theorem 4 (Hohenberg-Kohn). The mapping 



i H : [v] — ► P 1 H ([v]) — > -(f- 1 

pi 



(44) 



is injective for any H G "K p . 



Proof. It is sufficient to prove that [v] ^ [v'] =>- ih(M) ^h{[v'])- Let 
us suppose the contrary, that is [v] ^ [v f ] but ifl-(H) = ^h([ v ']) = P- Since, 
by definition, [v] ^ [v 1 ] implies W±{v) ^ W^(v'), it is possible to select 
a pair of eigenfunctions *f?i(v) G W^{v) and ^\{v') G W^{v') such that 
*&i(v ) G" Wi(v') and *&i(v') G" W/* (v ). Now, following the original arguments 
of Hohenberg and Kohn, we take into account that the quadratic functional 
associated with any Hermitian operator, reaches its absolute minimum at 
vectors belonging to the eigenspace of this operator corresponding to its 
lowest eigenvalue. We have two inequalities 

E?[v) < Tr {P?(v'){H + X{v))} = E?(v') + Tr {P?{v') [X{v) - X{v')}} 

= E^(v')+Tr{p(v-v')} 



and 



E?{v') < Tr {P?(v)(H + X{v'))} = E?(v) + Tr {P*{v) [X{v') - X(v)}}} 

= E? (v) + Tr{p(v' -v)} 

that contradict each other □ 

Let us point out the main differences between the original Hohenberg- 
Kohn statement [21 EHj and theorem 

(1) Theorem 4 states that the unique ensemble representable 1-density 
operator (not just density) corresponds to a certain class of 1-electron oper- 
ators, but not vice versa; 

(2) This theorem is valid for arbitrary choice of operator H that can be, 
in particular, arbitrary 1-electron, or even zero operator; 

(3) It is not presupposed that the mapping in is continuous (with respect 
to the quotient topology on Qh and the standard topology on Vn,p,i)- 

Note as well that the statement analogous to that in Theorem 0] holds 
true for the set "K q of q— electron operators (q < p) with obvious replacement 
of ^fC p_1 [^P 1 H (f)] in Eq. (j44j) by an ensemble representable density operator 
3jc p ~ q [j^Pi I { v )] °f order q. This generalization, however, is of no particular 
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interest because the structure of the set Vn,p,q of ensemble representable 
density operators of order q is unknown for q > 1. 

It is easy to show that in general case the image of the mapping (|44|) 
does not coincide with the set Vn, p ,i being only its proper subset. Indeed, let 
us suppose that the following auspicious conditions hold true (i) classes of 
1-electron operators from IKi can be continuously parametrized by free real 
parameters t = (t 1 , t 2 , . . . , t k ) G M. k : 

[v(t)] = {al n + v(t)} am (45) 

(/„ is 1-electron identity operator), and (ii) the mapping t — ► P^{[v(t)]) 
(see Eq.(jH])) is continuous (and, consequently, the mapping (jUj) is also con- 
tinuous). Even under such conditions the image of the mapping (|44|) is an 
open subset of the compact set Vn,p,i- It is pertinent to note as well that 
the mapping (J44j) is not necessarily continuous for arbitrary H. The simplest 
example corresponds to the case H = where the quotient set Q# is discrete. 

In the finite dimensional case the operator space IKi admits (in full anal- 
ogy with the set of representable 1-density operators, see Section II) at least 
two convenient realizations. In the first, standard, realization the MSO or- 
thonormal basis ip is supposed to be fixed and the operator space "K\ coincides 
with the set of all unitary transformed operators that are diagonal in the ba- 
sis under consideration. In the second realization a sum (disjoint union) of 
fibres isomorphic to the typical fibre W 1 is considered: 

% x = |J JCf (46) 

where "Kf is the set of 1-electron operators diagonal in the basis ip. It is clear 
that there exists natural surjective mapping !Ki — > !Ki. Of course, in the 
last realization the same 1-electron operator may appear on different fibres 
(for example, the 1-electron identity belongs to each fibre). This drawback 
is partially compensated by the relative simplicity of "Kx structure: it is 
homeomorphic to the Cartesian product the manifold N of all n-frames of 
1-electron space (base) and the typical fibre W 1 . The set "K\ is called a total 
space of a (trivial) fibre bundle with the manifold N as its base, JCf as a 
fibre over point ip, and MJ 1 as the typical fibre. 

In analogy with 1-electron operators, for each n-frame ip operator H 
has basis dependent realization with respect to the determinant basis 
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{ipi-t A ... A 4>i p }i L <...<i p - If if)' = 4>U then matrices of operators and 
are connected by 

= (k Uj H f (k uj (47) 

Of course, and H^i have the same eigenspaces being different represen- 
tations of the same operator. 

Since each fibre J-Cf is a subset of the set !Ki, the equivalence relation (f4*3j) 
induces the equivalence relation on !Kf for each n-frame if). The correspond- 
ing equivalence classes [v]^ constitute the set that can be considered as 
a fibre of the fibre bundle 

Qh=\J Q h (48) 

The replacement of the set Qh by the set Qh leads to the violation of the 
Hohenberg-Kohn theorem: the mapping 



(49) 



is obviously not injective. However, this theorem is valid for each fibre Jif, 
and even in more strong form: since 1-electron term Tr[pv] = Tr[d^(p)v] 
in average energy involves diagonal of representable 1-density operator, it is 
possible to prove that the mapping o i H is injective for any if). 

Thus, we can formulate the algebraic version of the Hohenberg-Kohn 
statement as 

Theorem 5. The mapping d o i H = {d^ o in}^^ is a fibrewise injective 
mapping of the fibre bundle Qh into the fibre bundle Vjv, p ,i- 

Corollary. Restriction of in on any fibre Q H is injective for any n-frame 

if,. 

Let us consider two in a certain sense boundary cases. 

For H = any fibre Qq is a discrete topological space with a finite number 
of elements (classes [v]^). Typical class [v]^ includes infinitely many diagonal 
operators 

n 
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with the following restrictions on their eigenvalues 

mak{nj} < /ip_fc + i = . . . = n p = . . . = Hi < min {//,} (51) 

3=1 3=1+1 



k,l 



n 



The ground state eigenspace of the corresponding operator X(v),v G [v 

is of the dimension ( l ~ p k hk )- By means of permutations of indices 1, 2, ... , 
applied to the typical inequality (J51|) it is easy to obtain n\/[(p — k)\(l — 
p + k)\(n — classes associated with the selected typical class. Elements 
of these classes are parametrized by n + p — I — k + 1 real parameters of 
which p — k parameters are free and the remainder ones obey inequality type 
restrictions. The total number of classes is 

^ / k=i Ky 7 i= P +i v 7 

For example, for iV = {1, 2, 3, 4} and p = 2 we have 

Six classes corresponding to single determinant ground states with typical 
inequality max{/^i,/x 2 } < mm{fi 3 , // 4 } and typical ground eigenspace C|12); 

Twelve classes corresponding to 2-dimensional ground eigenspace with 
typical inequality /ii < /i 2 = ^3 < fa and typical ground eigenspace C|12) © 
C|13>; 

Four classes corresponding to 3-dimensional ground eigenspace with 
typical inequality fa < fa = 7/3 = fa and typical ground eigenspace 
C|12) ©C|13) ©C|14); 

Four classes corresponding to 3-dimensional ground eigenspace with 
typical inequality fa = fa = fa < fa and typical ground eigenspace 
C|12) ©C|13> ©C|23); 

One class of the identity operator with typical inequality fa — fa = fa = 

fa. 

The image io(H]?j) of the typical class [v]^ is the density operator 



1 

P = ~ 

p 



,j=l j=p-k+l 



(53) 



Now let us consider the case when p-electron operator H defines the 
following structure of classes on the fibres Q^: 



n-l 



[v{t)f={aI n + Y,t 3 \^)m\ (54) 
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It is clear that this set of classes Qjj is homeomorphic to IR n_1 and each class is 
homeomorphic to R. If we suppose that the projector on the ground state of 
the operator H + X(v(t)) depends continuously on parameters £1, £2, • • • , tn-i 
then the image of the fibre is an open subset in Vtv, p ,i (or even in JV iPi i 
if the lowest eigenvalue of the parametric operator H + X(v(t)) is non- 
degenerate for all t). This means that with the class structure defined by 
Eq . (|54j) the image ij^(Q^) can not contain 1-density operators from the bor- 
der of Vjv,p,i (or Tat^i). All aforementioned arguments are applicable only 
if the ground state projector of operator H + X([v]^) is continuous (with 
respect to the quotient topology on Q%). And between these two bound- 
ary cases there exists a variety of 'intermediate' ones with more complicated 
structure of sets of classes. 

Parametrization of the set of classes defined by Eq. (j54|) is by no means 
unique. Sometimes it is convenient, for example, to impose the additional 
requirement of orthogonality of the parametric representative of class [v (t)]^ 
to the identity operator (say, with respect to the trace scalar product). It 
may be easily done by writing class [v (t)}^ in the form 

[v{t)f =\al n + j^t 3 \^)(^\\ (55) 
^ i =1 > am 

with the additional condition ^2 tj = 0. Concrete examples of different 

jeN 

parametrizations will be given later. 

Now let us turn to more detailed analysis of the 'diagonal' version of 
the Hohenberg-Kohn statement and consider the mapping o i H : — > 
VjVpi for some fixed n-frame ip. It will be supposed that matrix of p- 
electron operator H has non-degenerate ground state. Parametric matrix 
H^ + X (v^(t)) either has non-degenerate ground state for all finite values 
of parameters, or there are points in the parameter space where the ground 
state is degenerate, or quasi-degenerate. In strictly degenerate case, due to 
Theorem 5, the mapping in should have discontinuity when going from pure 
to averaged 1-density diagonal. 

It should be specially emphasized that the structure of classes from 
depends on the properties of representable diagonal d^{p) of 1-density oper- 
ator p corresponding to the ground state of operator H. In particular, this 
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structure is different for d^(p) G V% p i and d^{p) G cA/^ v Indeed, let 

fo Hie I 

"Hi Hie J (56) 
l p 

where I (1 J = 0. This means (see Appendix A) that <i,/,(p) belongs to the 
(/, J)-face of the polyhedron V% p i of dimension n — \I\ — \ J\ — 1. Ground 
state wave function of matrix corresponds of 1-density operator with 
diagonal elements satisfying conditions (JSTIj) if it is orthogonal to subspace 
spanned by p-electron basis determinants \R) such that either R n / 7^ or 
R fl J = 0, and arbitrary matrix of the form 

min{|/|,p} (k) (p—k) 

H^a,b)=H^ + aI n + J2 b KyL \K U L)(K U L\ (57) 

k=0 Kcl LcN\(IUJ) 

has as its eigenvector. Here is 1-electron identity matrix. To describe 
the structure of the set Q^, it is necessary to find out under what restrictions 
on the coefficients bx,L 

(i) the parametric term on the right hand side of Eq. fJSTj) is 1-electron 
matrix of the form of Eq. (|41j) , 

and 

(ii) has as its lowest eigenvector. 

Let us consider an example formally corresponding to 2-electron case: 
N — {1,2,3,4} andp = 2. 

For a fixed 4-frame ip = (ipi, ip2, ^3? ^4) general two-electron wave function 
is of the form 

^ = aipi Aip2 + Hi A ip 3 + cipi At/> 4 + dip 2 A ip 3 + eip 2 A ^4 + fip3 A ^4 (58) 
Diagonal of 1-density operator p corresponding to this wave function is 

d^p) = l -{a 2 + b 2 + c\ a 2 + d 2 + e 2 , b 2 + d 2 + f 2 , c 2 + e 2 + f 2 ) (59) 

We analyze three cases: (1) d^(p) belongs to the interior part of the polyhe- 
dron y%2i'i (2) d^p(p) belongs to a 2-face of this polyhedron, and (3) d^{p) 
belongs to an edge of the polyhedron V^ 2 x . 
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With the aid of MATHEMATICA 4.1 package [T5] three 6 x 6-matrices 
1 , and with predefined spectrum and eigenfunctions, corre- 
sponding to these three cases, were generated (see Appendix C). 
General form of 1-electron operators belonging to the fibre "Kf is 



H (i) 



'(t) 










h 
o 
o 



o 
o 

h 










(60) 



and the corresponding 2-electron operator may be written as (see Eq. (l41|) 



X (v*(t)) 



ft x + 1 2 





V o 





h + t 3 
o 
o 
o 


















t 2 + h 











t 2 + U 




o \ 






f 3 + Uj 



(61) 



Case 1. d f (p) = (0.404545,0.313636,0.163636,0.118182) G V 



AT,2,1 - 



The lowest eigenfunction ty^ of the matrix H^' involves all 6 determi- 
nants and all coefficient bjc,L i n Eq. (j57|) should be taken equal to zero. This 
is the most general case and classes of 1-electron operators are of the form 



ah + 



(x 

y 

z 

\0 









\ 



(62) 



-x — y — z J 



aelt 



where x = 3t\ — 1 2 — 1% — £4, y = 
Each class is homeomorphic to 



-ti + 3t 2 — £3 — 1 4 , and z = —ti —t 2 + 3t^ — 1$ 



and the set of classes is homeomorphic 



to 



Note that in this concrete case the set of classes carries natural vector 



space structure. 

An example of behavior of diagonal elements of 1-density operators along 
one of basis directions in the parameter space is displayed on Figsffl It is 
seen that these elements are smooth functions along the selected direction 
and that the mapping z — > (pn(0, 0, z),p 22 (0, 0, z), p 33 (0, 0, z), pu(0, 0, z)) is 
injective. It may be expected that each representable diagonal corresponds 
to the unique set of parameters, and, consequently, to the unique class of 
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Figure 1: Diagonal matrix elements of parametric 1-density matrix p(x, y, z): 
functions z — > pu(0, 0, z) 

1-electron operators. It is also clear that for any finite values of parame- 
ters the corresponding diagonals of 1-density operators belong to the interior 
part of the polyhedron 2 x in full accordance with the standard topological 
arguments. It seems to be most likely that for the operator under considera- 
tion the (topological) closure of the image of the set with respect to the 
mapping o i H coincides with the polyhedron V^ 21 . 

Case 2. d^(p) = (0.333333,0.333333,0.333333,0) G F {4}j0 where F {4}j0 = 
Conv (t)2j.i(12), ^24.1(13), U2j.i(23)) is 2-face of the polyhedron V^ 21 (see Ap- 
pendix A). 

(2) (2) 

The lowest eigenvector ty\ of matrix H\ is orthogonal to basis deter- 
minants 1 14), 1 24), and |34). However, the parametric operator (see Eq. (|57j) ) 

4 2) (a, 6) = Hf + \aX{I A ) + 6X(J 4 ) (63) 
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where 



J A 



/o 











(2) 

has ty\ as its lowest eigenvector only for bo < b < oo where bo s 
(see FigO). As a result, at the point bo the class structure undergoes essential 



(64) 



-5.5484 




Figure 2: Two lowest eigenvalues of the parametric operator H^'(a,b) (see 
Eq.(jnHl)) as functions of parameter b 
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change. We have 



ah + b Ja + 








0\ 1 





a/4 + 



(x 

y 

-x-y 

\0 0/ J 

/x oV 

y 

-x-y-b 

\0 b) J 



> if 6 > 6 

(a,b)eKx (60,00) ^g[^ 

if 6 < 60 



Case 3. d^(p) = (0.5,0.470588,0.0294118,0) e %h{i} where F {4}i{1} = 
Com; (t> 2 ji(12), 1)24.1(13)) is 1-face (edge) of the polyhedron V]y- 21 (see Ap- 
pendix A). 

The lowest eigenvector \E^ 3 ' ) of matrix if^ 3 ' ) is orthogonal to basis deter- 
minants 1 14), 1 23), 1 24), and |34). The auxiliary parametric operator 



H®(a, b, c) = H ( ;> + ^-aX{h) + bX{J 4 ) + cX(K A ) 



r(3) 



where 



(\ 0\ 
0-100 

0-10 

\0 0/ 



(66) 



(67) 



(3) 

has ^ as its eigenfunction. To determine the relevant class structure, it is 
necessary to find out when the parametric operator (J66)) has tyf^ as its lowest 
eigenvector. From FigEI it is seen that the set of admissible parameters is 
a certain subset S of the plane E = 10 in 3-dimensional space of triples 
(6, c, E). Within this subset the class structure is 



[v{x)\ 



clIa + bJA + cKa + 



/0 0\ 
0x00 
-x 
\0 0/ 



(a,6,c)eKxS 

After passing through the border of the set 5* the class structure changes. 
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Figure 3: The lowest eigenvalue of the parametric operator (a,b,c) (see 
Eq.fjnni)) as a function of parameters b and c 

Thus, sets of classes of 1-electron operators associated with some p- 
electron operator may be of a rather complicated nature. It may be 
a union of cells of different dimensions as it was in the last two examples, 
and use of general parametric classes of the type of Eq. (j55)l for both interior 
and border points of the polyhedron V% p 1 will lead to violation of Theorem 
5. Probably the only regular case corresponds to the situation when repre- 
sentable diagonal of 1-density operator, corresponding to the ground state 
of operator H, belongs to the interior part the Coleman's polyhedrons , 
for any n-frame ip (see Proposition 2). In this case the following condition 
may be fulfilled: 



for any n-frame ip. If the equalities ()69|) hold true, then the inverse mappings 



are correctly defined. 

Till now we studied in a certain sense local task of parametrization of 
represent able diagonals by classes of diagonal 1-electron operators from a 
given fibre "Kf . The next very important step is to assemble local results. 




(69) 



(70) 
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Since the fibre bundles under consideration are trivial (that is not twisted, 
in contrast to, say, the famous Mobius band), this task, in theory, is not 
complicated. Choosing any convenient (local) parametrization of the base 
N (dimN = n(n — l)/2 if the ground number field is K), we come to 1- 
density diagonals depending on two independent sets of parameters. The 
first set parametrizes some open subset of the base N and the second set 
parametrizes classes of diagonal 1-electron operators. 

Now let us consider the case when for each n-frame if) the structure of 
classes of 1-electron operators from Jif is given by Eq.()55|) without explicit 
reference to any concrete p-electron operator H. Each fibre (subscript 
H is omitted because of the aforementioned reason) is (n — l)-dimensional 
vector space over K with respect to the operations 

[v(t)]^ + {v{t')f = [v{t + t')f, a[v{t)^ = [v{at)f, aeR (71) 

and the corresponding fibre bundle Q together with the vector space structure 
on each fibre is the so-called vector bundle JTj . The following problem seems 
to be of primary interest for the DFT theory: 

Under what conditions the inverse of a fibrewise bijective mapping 

s ■ U - U Q " ( 72 ) 

is the Hohenberg-Kohn mapping d o i H for some p-electron operator HI 
It is easy to present the following very important necessary condition that 

should be imposed on a fibrewise mappings of the type of Eq. (j72Jl : 

For any fixed n-frame if) and any 1-electron unitary transformation U 

Sipu {d^u {U ] p^U)) = W u (73) 
where is some 1-density operator belonging to the interior part of the set 

In concluding this section let us discuss in general form the universal 
functionals of the DFT theory (see [El EI and references therein). Relevant 
definitions and theorems from set-theoretical topology may be found in [TT] . 
Let us suppose that 

7 : P(^, p ) - 7 (74) 

is a continuous surjective mapping of the (projective) space of states on a 
Hausdorff topological space T. Being the image of the compact set P(JF /Vp ) 
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Figure 4: Branches of the universal functional on a fibre P]y 2 i corresponding 
to the natural MSO 4-frame of operator HjP ground state 

with respect to continuous mapping, T is also compact. Let i? 7 be the equiv- 
alence relation associated with the mapping 7: 

Rr(x,y)^-y(x) = j(y) (75) 

The quotient ¥{Tn,p)/ R-i is necessarily Hausdorff and compact space (in quo- 
tient topology) and, consequently, the mapping 7 admits representation in 
the form 7 = gov where n is the canonical projection P(jF A r iP ) — > P( J~N,p) /Ry 
and g is the homeomorphism P( J-N,p) /Rj - > 7. The average energy 
Eh = (9\H\ty) is a smooth mapping P(jF/v iP ) — > R. Of course, this mapping 
is not compatible with the equivalence relation i? 7 . But it is possible to 
introduce a new mapping E H as 

E H :-y-\r)^ inf E H (x), re7 (76) 

x£-y 1 (r) 

that is correctly defined due to the Weierstrass theorem [T2j. The universal 
functional may be written as F = En o g^ 1 ; 7 — ► M. It is continuous if the 
mapping E H is continuous. 

If 7 = ^c p ~ 1 then 7 = ^Pn,p,i- The domain of the universal functional 

F : ?jv, p ,i — > K, even if the mapping is continuous, not necessarily admit 
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global parametrization (see Appendix B). Therefore, formal manipulations 
involving derivatives and variations of universal functionals, require certain 
care. 

Let us suppose that the set JV, P ,i is realized as a fibre bundle (|T9"j) and 
consider simple example: N = {1,2,3,4} and p = 2 (see Appendix B). In 
this case there are three branches of universal functional associated with 
three pairs of disjoint subsets: F R}N \ R (X; ip). If 2-electron operator H^p (see 
Appendix C) is chosen, then the ground state branch corresponds to the pair 
(12,34) (it is sufficient to find the ground state 1-density operator and then 
perform transformation to the natural MSO basis to see that in this basis 
ifW = 0.944996|12) - 0.327082|34)). Dependence of 3 branches of the the 
universal functional on the parameter r (see Appendix B, Eq.(B.2)) on the 
fibre Pjy-2 ij corresponding to the ground state natural 4-frame tp, is displayed 
on Fig. I3J In the case under consideration these branches are just segments 
of three ellipses. 

V. Conclusion 

Any finite dimensional model of quantum mechanics is in a certain sense 
incomplete, because it does not include many important relations like com- 
mutation relations of canonically conjugate pairs of observables which exist 
only in infinite dimension. On the other hand, in the majority of cases molec- 
ular calculations are performed in finite basis sets, and it seems reasonable 
to try to find out what specific finite dimensional features are inherent in 
algebraic version of quantum chemistry methods. Thorough analysis of the 
algebraic version of the Hohenberg-Kohn theorem for arbitrary selected many 
electron operator H shows that 

- structure of classes of 1-electron operators that appear in the Hohenberg- 
Kohn theorem may be very complicated and may be determined not only 
by operator H ground state but also by spectrum of a certain auxiliary 
parametric many electron operator; this structure strongly depends on the 
properties of the operator H ground state 1-density operator p and is different 
for p belonging to the border of the Coleman's set and to its interior part; 

- if p belongs to the interior part of the Coleman's set (all MSO's are nec- 
essarily active) then there exist many electron operators such that the corre- 
sponding classes of 1-electron operators are of general form {al n + v(i)} a&R ; 
in this case the Hohenberg-Kohn mapping may parametrize the interior part 
of the Coleman's set but its image can not contain, say, vertices of this set 
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(HF 1-density operators); 

- the Hohenberg-Kohn mapping is not necessarily continuous for arbitrary 
operator H; 

- in finite dimensional case fundamental role is played not by densities 
but by representable diagonals of 1-density operators; 

- the set of 1-density operators representable by pure states, in general, 
does not admit global parametrization and there may exist several branches 
of DFT universal functionals. 

General theory developed in this work does not impose any special re- 
strictions on many electron operator H. However, in all examples considered 
it was assumed that the ground state of this operator is non-degenerate. De- 
generate case may introduce additional complications in concrete structure 
of classes of 1-electron operators and requires separate investigation. 
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Appendix A: Combinatorial Structure of Polyhedron Vn, p ,i 

Polyhedron Vjv, P ,i is situated in the (affine) hyperplane 

H n = {x e R n : (a\x) = 1} (Al) 

of the vector space M. n that has the vector 

n 

i=i 

as its normal. Here {e^^N is the canonical basis of lR n and N = {1,2, ... ,n}. 

It is convenient to recast the Coleman's system (10) in a form commonly 
accepted in theory of polyhedral sets: 

f (e,|A) > 0, i G N 

< (a-pei\\) > 0, ieN (A3) 
U|A) = 1 

To describe faces of the polyhedron Vn,p,i h is necessary to analyze 'mixed' 
families {{ej} ig /, {o — pej}j € j, a} with I D J = 0. Two obvious restrictions 
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should be imposed on disjoint subsets / and J: < \I\ < n — p and < 
| J | < p. Indeed, if \I\ > n — p then the number of non- vanishing components 
of vector A G Vw,p,i is less than p and consequently, (ci| A) < 1. If |J| > p 
then the number of components of vector A G Vjv, p ,i equal to ^ is greater 
than p and, consequently, (a|A) > 1. Direct calculation shows that a family 
{{ei}iei, — P e j}jeJ, &} is f ree if an d only if Id J = and |/U J\ < n. Thus, 
faces of the polyhedron Vn, p ,i are the sets F I:J of solutions of the systems: 

(e,|A) = 0, tel 
{a-pe^X) = 0, j G J 

(a|A) = 1 (A4) 
(e;|A) > 0, i G iV\J 
(o-pej|A) > 0, j G iV\J 

Let /, J be subsets of the index set N, such that Id J = 0, \IU J| < n — 1, 
< |/| < n - p, < | J\ < p, and let K = N\(I U J). Let us put \I U J| = r 
and | J| = p — s. Admissible values of r are 1, 2, . . . , n — 1. For fixed value of r 
the admissible values of |J| are min{r, p}, min{r, p} — 1, . . . , max{0, r — n+p}. 
If A G then the following equality holds true 



1^1 
p 



+ $> fc = l (A5) 



k&K 

This equality can be recast in a more convenient form: 

s 



5> 

keK 



(A6) 



where s = max{p — r, 0}, max{p — r, 0} + 1, . . . , min{p, n — r}. There are 
three cases to be analyzed. 

(1) s = max{p — r, 0} = 0. 

In this case Eq.(A.6) has only trivial solution A& = 0, k G K. The 
corresponding face F^j is just the vertex J) (0-face). 

(2) s ^ and |X| > s + 1. 

The set of solutions of Eq.(A.6) is non-degenerate Coleman's polyhedron 
|Vk",s,i of the dimension \K\ — 1 and, consequently, (\K\ — l)-face of the 
polyhedron \/n, v ,i- 

(3) s 7^ and \K\ = s = n — r. 
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In this case we have degenerate Coleman's polyhedron ~\/k,s,i — 
{^v s ii(K)} of dimension (vertex) and, consequently, 0-face f p ji(J U K) 
of the polyhedron Vn, p ,i- 

Now we can calculate the number of faces of a given dimension n — r — 1 
of the polyhedron Vn, p ,i- 

((;) if r = n-l 

f -, — / min{p,n— i — 1} ( A 7\ 

^ s=max{p-r,l} 

For example, for n — 12 and p = 5 the /-vector (see, e.g. [16j) is / = 
(1, 792, 13860, 46200, 76230, 77616, 52668, 24552, 7920, 1760, 264, 24, 1) where 
we have added two coordinates /_i = 1 and / n _i = 1 corresponding to 
improper faces and V/v,5,i- It is easy to check that this /-vector satisfies 
the classic Euler identity: 

n-1 

£(-!)*/* = (A.8) 

k=0 

Summing up, we can state that each (n — r — l)-face of dimension greater 
than of the polyhedron Vjv,p,i is defined by a pair of disjoint subsets (J, J) 
with | J U J\ = r excluding the cases \J\ = p and | J\ = r — n + p. A vertex 
Vpii(R) belongs to this face if and only if J C R C N\I. The number of 
vertices belonging to such a face is equal to • Coleman's polyhedrons 

form a chain 

I a j = VjV, n ,l C Vjv,„-l,l C . . . Vat.2,1 C Vat,!,! (A9) 

where Vjv,i,i and V^n-i.i are standard and non-standard simplexes of the 
vector space R n , respectively. 

Appendix B: An Example of Geometric Description of the Set 
of 1-Density Operators Representable by Pure States 



Let us consider general two-electron wave function of the form of Eq. 
The corresponding pure state belongs to the projective space P^at^)- 

We confine our consideration to the case when the ground number field is R, 
and the projective space can be realized as a quotient of the 5-dimensional 
unit sphere S 5 formed by pairs of diametrically opposite points. 
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Matrix of 1-density operator with respect to fixed basis ip obtained by 
contraction of is 



(a 2 + b 2 + c 2 bd + ce —ad + cf —ae — bf\ 

1 



P= 2 



bd + ce a 2 + d 2 + e 2 ab + ef ac — df 
-ad + cf ab + ef b 2 + d 2 + f 2 be + de 
\ -ae -bf ac - df be + de c 2 + e 2 + f 2 J 



(B.l) 



Let us suppose that 1-density matrix (B.l) is diagonal in the basis under con- 
sideration. We have the system of 6 polynomial equations with respect to 
variables a,b,c,d,e, f that can be easily solved with the aid of MATHEMAT- 
ICA built-in Reduce function ^S] • For each pair of 2-electron determinants 
with disjoint index sets there is a solution 

p(R,N\R) = rv 2ll (R) + (1 - r)v 2il (N\R) (B.2) 

where R is 2-element subset of the index set N. 

Using iteration formula (|34|) it is possible to construct infinitely many 
(for r ^ 0,1) 2-electron determinant ensembles contracted in diagonal (B.2). 
The corresponding pure states (see Proposition 1), however, give diagonal 
1-density operators (over IR) only if 

V ± = t?\R)±(1-t)*\N\R) (B.3) 

Thus, in the realization of the set ?V,2,i as a fibre bundle with the mani- 
fold of 4-frames as its base the fibre 2 1 over point if) consists of 1-density 
operators of the type of Eq.(B.2). The number of different types of such 
operators is equal to 3. In geometric terms the fibre P% 2 i ^ s a um° n of 3 
closed line segments having exactly one central point (|, t, |, |) in common. 
Representable by pure 2-electron states 1-density operators, diagonal with 
respect to the basis if), necessarily have the form of Eq.(B.2), and, conse- 
quently, dimP^ 21 = 1. CPjv,2,i may be realized as a union of three (trivial) 
fibre bundles 

3 

?W = U ^2,1 IP A) 

i=l 

where each component (pjy 2 i i s homeomorphic to the Cartesian product N x 
[0,1]. 
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Appendix C: Construction of Parametric Diagonal Elements of 
Representable 1-Density Matrices with MATHEMATICA Package 



Matrix H$ } was constructed using the following code: 

In[l] :=<< Linear Algebra" Orthogonalization" 
In[2] := {ul, u2, u3, u4, u5, u6} = 

GramSchmidt[{{0.8, 0.4, 0.3, 0.2, 0.1, 0.4}, {1, 0, 1, 0, 0, 0} 
{0,1,0,1,0,0}, {0,0,1,0,1,0}, {1,0,0,0, 0,1}, 
{-1,1,1,-1,1,1}}]; 
In[3] := ul = Outer [Times, ul, ul]; u2 = Outer [Times, u2, u2]; 

u3 = Outer [Times, u3, u3]; u4 = Outer [Times, u4, u4]; 
u5 = Outer [Times, u5, u5]; u6 = Outer [Times, u6, u6]; 
H = -10 * ul + -5 * u2 - 4 * u3 - 3 * u4 - 2 * u5 - u6; 



Diagonal elements of 1-density matrix as functions of free parameters x, y, z 
(see Eq.(jnni)) were defined as 

Jn[4] := pn[x_, y_, z_] := Block[{m, v, h, ev, vec, ind, a, b, c, d, e, f , t}, 
H = -10 * ul + -5 * u2 - 4 * u3 - 3 * u4 - 2 * u5 - u6; 



fx + y 

x + z 












y + z 



















Xv 









-y-z 









— x — z 



\ 







x-yj 
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h = H + Xv; {ev, vec} = Eigensystem[h]; 
ind = Ordering[ev]; k = ind[[l]]; 

a = vec[[k, l]];b = vec[[k, 2]];c = vec[[k, 3]];d = vec[[k,4]]; 
e = vec[[k, 5]]; f = vec[[k, 6]]; t = 0.5 * (a 2 + b 2 + c 2 )]; 
Plot[p n [0,0,z],{z, -50,50}]; 



Matrices H^> and H^> were constructed in the analogous manner with the 
initial linearly independent vectors 

(^,0,-^,0,0), (0,1, 1,0, 0,0), (0,1, 0,1, 0,0), (0,0, 1,0, -1,0), 
(1,0,0,0,0,1), (-1,1, 1,-1,1,1), 

and 

(|, 1 0, 0, 0, 0), (0, 1, 1, 0, 0, 0), (0, 1, 0, 1, 0, 0), (0, 0, 1, 0, -1, 0), 
(1,0,0,0,0,1), (-1,1, 1,-1,1,1), 

respectively 
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